Dynamic response of trapped ultracold bosons on optical lattices 
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We study the dynamic response of ultracold bosons trapped in one- dimensional optical lattices 
using Quantum Monte Carlo simulations of the boson Hubbard model with a confining potential. 
The dynamic structure factor reveals the inhomogeneous nature of the low temperature state, which 
contains coexisting Mott insulator and superfluid regions. We present new evidence for local quan- 
tum criticality and shed new light on the experimental excitation spectrum of 87 Rb atoms confined 
in one dimension. 
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PACS numbers: 03.75Hh,03.75.Lm,05.30. Jp 

The advent of Bose-Einstein condensation in cold 
atomic gases opened new fields of research merging 
atomic physics and quantum optics with condensed mat- 
ter physics. Scattering experiments common in solid 
state physics and liquid helium have been adapted to 
extract information on excitation dynamics of atomic 
condensates resulting in the technique of Bragg spec- 
troscopy 0] . Initial experiments on Bose-condensed gases 
could be well described by mean-field Gross-Pitaevskii 
approaches, valid for weak interactions. However, recent 
experiments have reached the regime of strongly inter- 
acting atoms, often by confining the atoms in optical lat- 
tices, leading, for instance, to the destruction of a phase- 
coherent superfluid (SF) phase and the establishment of 
Mott insulator (MI) phase with squeezed number fluctua- 
tions (2,13, El- The latter state is considered as promising 
for use as a large quantum register for quantum informa- 
tion purposes |5|. 

Ultracold atoms in optical lattices are an almost ideal 
realisation of the boson Hubbard model [f| , that has pre- 
viously been studied intensely (with different motiva- 
tion) 0, 13- Much of this previous work focussed on 
the equilibrium phase diagram or on the conductivity 
in the presence of disorder. However, as in other con- 
densed matter systems, a rich source of information on 
the excitation dynamics is the dynamic structure factor, 
the double Fourier transform of density-density correla- 
tions, which can be measured using Bragg spectroscopy 
HI EJ lioj - It is therefore of great current interest to 
study the dynamic structure factor of the boson Hub- 
bard model. Because of loss of translational invariancc 
due to the presence of the confining trap, results for the 
uniform system cannot be directly applied to the con- 
fined system. This affects both the static response, the 
system no longer has a globally incompressible phase [Tl). 
and the dynamics, since the wave vector is no longer a 
good quantum number. 

In this paper, we use Quantum Monte Carlo to study 



the ground state of the one-dimensional boson-Hubbard 
model in a confining harmonic potential: 
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where L is the number of sites and Xi is the coordinate of 
the ith site. The hopping parameter, t, sets the energy 
scale, rij = a\ai is the number operator, [oj,at] = 
are bosonic creation and destruction operators. Vr sets 
the confining trap curvature, while the contact interac- 
tion is given by U. We use the World Line algorithm 
for our simulations. The static properties and state di- 
agram of this model were studied in (lll | where it was 
shown, among other things, that the presence of the trap 
and the destruction of translation invariance eliminate 
the quantum phase transitions between the SF and MI 
phases and replace it by a co-existence of phases. Our 
main results are: (i) due to co-existence of MI and SF re- 
gions in the confined system, the dynamic structure fac- 
tor exhibits two excitation branches, one gapped and one 
with a linear, phonon-like, dispersion; (ii) these branches 
are shown clearly to be connected to the spatial regions 
containing the MI and SF phases and (iii) by using an 
appropriate measure of local excitations, we provide new 
evidence for local quantum criticality. 

In the ground state (T = 0) the dynamic structure 
factor is given by 
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where Lt is the number of time slices, E m the energy of 
state |m), and (n(r, 0)n(r', t)) the space and time sepa- 
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FIG. 1: (a,c) S(k,u>) vs lo and (b,d) fi(fe) vs k for k — 2im/L 
with 1 < n < 9. The system has commensurate filling: L = 
20, N b = 20,/3 = 10 with (7 = At (a,b) and [7 = It (c,d). 



rated density-density correlation function. With our nor- 
malization we have J dcoS(k,Lo) — NbS(k) where S(k) is 
the static structure factor and Nj, the number of bosons. 
Equation shows that, for a given k, the peaks in 
S(k,uj) correspond to the excitation energies. We have 
performed all our simulations at /3 = T _1 = 10 and veri- 
fied that it is sufficient to study the ground state. In the 
World Line algorithm, it is a simple matter to measure 
(n(r, 0)n(r', r)) where r is the imaginary time separation. 
The dynamic structure factor is then given by Eq. J5J, 
with t replaced by —it which then requires the perfor- 
mance of a Laplace, instead of a Fourier, transform. This 
is done with the help of the stochastic method (SM) . 

Figure n shows S{k,Lo) for the uniform system in the 
MI (a) and SF (c) phases. The critical value in one 
dimension is U cr i t — 2.33i 0]. For better visibil- 
ity, our S(k,Lo) curves in Figs. IT131 have been normal- 
ized J2uj S(k,cu) = 1 and do not show the true spectral 
weights. For the SF, it is clear that the excitation ener- 
gies for small k are linear in k, i.e. phononic, leading to 
a stable SF phase (according to the Landau criterion). 
On the other hand, the MI case clearly shows a gap at 
lo « it, where the spectral weight in Fig. ^a) starts to 
be appreciable. This is in excellent agreement with the 
value of the gap of 3i measured in [8J from the static 
compressibility. Note that the first peak is at lo » 4.5t, 
which is larger than the gap, and that the peaks for the 
MI are rather wide. This is due to the presence of sev- 
eral closely spaced peaks which we cannot resolve. Such 
closely spaced peaks are seen in the exact diagonalization 
of smaller systems where, for the MI, the first in the 
lowest group of peaks does indeed correspond to the gap. 



One can study the excitation spectrum using the dis- 
persion relation fi(fe), 



n(jfe) 



J duj lo S(k, u>) 
J du> S(k, ui) ' 



(4) 



which can be expressed in terms of purely static quanti- 
ties using the /-sum rule 



do; lo S(k, lo) = N(,Ek, 



(5) 



E k = -t (cos(27rfc/L) - 1) (0| (4^+1 + 4+iOi) 1°) 



(6) 

The dispersion relation is then given by the Feynman 
result, 



E k 

S(k)- 



(7) 



Figure n shows fi(fe) obtained using the dynamic struc- 
ture factor, Eq. J3J (circles), and also using only static 
quantities, Eqs. © and J7J) (pluses) in the MI (b) and SF 
(d) phases. These two calculations agree very well and 
give for the SF phase linear dispersion for small k and 
a gap for the MI. While £l(fc) shows clearly the presence 
of a gap for MI, it does not give its value since it is an 
integral over all peaks that may be present. In addition, 
we show the positions of the peaks in both cases (down 
triangles) which agree extremely well with Q (k) when the 
peaks are very sharp. This offers a consistency and pre- 
cision check on our methods. Our results for fi(fc) agree 
with the small lattice exact diagonalization results [l0( . 

We now turn to the confined system where we always 
take a system size L = 100 and Vt — 0.008. These values 
and the fillings we use are comparable to the experimen- 
tal ones. Figure shows the local density profiles for 
a system in a pure SF phase, (a), and with coexistence 
of SF and MI, (b). Figure |3| shows the corresponding 
S(k,oj). For the pure SF case, Fig. G^a), S(k,uj) is sim- 
ilar to the SF case in Fig. QIc): no gap is visible and 
the excitation energies go to zero smoothly as k — > 0. 
Figure E{b) shows the same for iV& = 64 where there is 
coexistence of MI and SF phases. There is a marked dif- 
ference in the behavior of S(k,u>) compared to the pure 
SF case: two branches of excitation peaks are clearly vis- 
ible in the form of two ridges. The low energy branch is 
seen to be very similar to S(k, uj) in the SF case, Fig.^a), 
saturating at about u) « M, and going smoothly to zero. 
The second family of peaks saturates at a higher energy, 
lo w 8t and seems to make its first appearance at about 
k w 11 in units of 2-k / L with anw« 3t. 

To interpret these two excitation branches, we ex- 
amine S(x, lo), the imaginary time Laplace transform 
of {n(x, 0)n(x, r)), the same site but time-separated 
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FIG. 2: Local density profile, p(x), for the confined system 
with L = 100, V T = 0.008, U = At: (a) N b = 25 (pure SF) 
and (b) N b = 64 (co-existence SF-MI). 
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FIG. 4: S(x,u>) for the system in Fig. E^b) and shows the 
origin of gapped and ungapped excitation branches are the 
SF and MI regions. 
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FIG. 3: S(fc,w) vs lu (k = 2-rvn/L) for the system in Fig. 
N b = 25 (a) and N b = 64 (b). The two ridges in (b) reflect 
the SF-MI co-existence, Fig.^b). 



density-density correlation function. This is shown in 
Fig. 01 which clearly identifies excitation branches with 
position in the trap. S(x,uj) includes contributions from 
all k values. The very low energy phonon excitations 
arise from the SF regions in the outer wings and the cen- 
ter of the system. The middle excitations, uj ~ 2t — 3t 
appear prominently in the squeezed SF region between 
the MI and the outer edges of the system. The value of u 
at the peak corresponds well to the saturation value, at 



high k, of the SF branch of S(k,u>) in Fig. CD This effect 
was observed in all our simulations where a MI phase is 
present. In these outer regions the superfluid is squeezed 
between the MI and the rapidly increasing confining po- 
tential which appears to enhance the spectral weight of 
the high uj excitations and make them so prominent. Fi- 
nally, the peaks in the two MI regions in Fig. 0] agree 
well with the saturation value of the higher uj branch 
in Fig. 01b) . This confirms our interpretation that this 
branch corresponds to excitations in the Mott phase. As 
further evidence, we mention that the S(x, uj) for the uni- 
form system in the MI phase at U — At (Fig. ^a)) has 
its peak atw« 8t, like the Mott regions in Fig.|3Jb). 

The interpretation of the higher uj excitation branch in 
S(k, uj) as being due to the MI regions (Fig. 0} and the 
absence of a gap in Fig.[3]show that even when MI regions 
are present, the system is gapless! Figure 0] shows that 
local excitation energies in the MI regions correspond to 
those in MI phase of the uniform system at the same U. 

Another very striking feature in Fig. 21 is the behavior 
of S(x, uj) at the boundary between SF and MI sites: The 
peak in S(x,uj) is very broad and extends from uj = to 
the largest values examined (uj = 20t) which indicates a 
diverging correlation length in the time direction, £t(ir). 
A diverging £ t means that (n(x, t)ti(x, 0)) decays as a 
power law in r and characterizes a quantum phase tran- 
sition. In this case, since it happens only locally, the 
divergence of £t (x) is renewed evidence of the local quan- 
tum criticality already discussed in [ill [LU . 

To measure the total response of the system, summed 
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over all momenta, we use 

S(uO=£S(fc,w). (8) 

k 

This quantity is shown in Fig. [S] for the system shown in 
Fig. Efb) (same as Fig. EJ) and also for N b = 50, U = It 
which has pure MI in the center surrounded by SF at 
the edges. A two peak structure is seen where the first 
peak appears not to depend on U whereas the second 
peak does. The first peak, therefore, is due to the SF 
component of the system and has an oj which corresponds 
very well to the saturation value of the SF branches of 
S(k,oj) in Figs. 0Ja,b). O n t ne other hand, the higher- 
peaks depend on the value of U and both correspond to 
the saturation values of the MI branches in S(k, oj) (see 
Fig. PJV)). Both peaks come from the high k values of 
the SF and MI branches because those values have the 
larger spectral weights which therefore dominate when 
the sum over k is taken. Recall that the S(k,oj) in the 
figures have been rescaled for better visibility and do not 
show the true spectral weights, but Fig. |S] does. 

Our two-peak structure, Fig. resembles the ex- 
citation spectrum in Fig. 2(a) in {Sj for 87 Rb on one- 
dimensional lattices. Both quantities measure the total 
response of the system to a given frequency, but they are 
not the same and the connection between them is not 
clear. Although S(k, oj) does detect the presence of MI 
regions, it also shows there is no overall gap but that lo- 
cal excitation energies do correspond to gap energies in 
the uniform system at the same U. It is possible that 
the SF peak we detect is missed by reference because 
the amplitude with which they shake their system is very 
large (20% of the optical lattice potential). 

We have presented the first Quantum Monte Carlo 
study of the dynamical structure factor, S(k,oj), for the 
Hubbard model on uniform and confined optical lattices. 
The uniform system exhibits linear dispersion (phonons) 
for the SF and a clear excitation gap in the MI in agree- 
ment with the directly computed gap || . In the absence 
of MI regions in the confined system, the S(k, oj) behavior 
is very similar to the uniform case. 

However, when MI regions are present (Figs.^b) and 
0}, we find that although S(k,oj) shows the presence of 
such regions, it also clearly shows the absence of a gap: 
the system is globally compressible and gapless . The 
presence of the MI regions is deduced from the presence 
of two excitation branches in S(k,oj) or from the double 
peak structure of S(oj) (Fig. [SJ. We note, as discussed 
above, that the second peak in Fig. El corresponds to MI 
regions. In addition, we note that S{k,oj) and S(oj) give 
information about local excitations in the trap. For ex- 
ample, by comparing S(k, oj) with S(x, oj) (Figs.[3{b) and 
0} we have found that the first peak in S(oj) represents 
excitations in the SF regions at the edges of the system 
while the second peak represents the excitations in the MI 



regions. The value of the local excitation energy is con- 
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FIG. 5: S(oj) for N b = 64, U = At and for N b = 50, U = It 
which has only a MI in the middle. 



sistcnt with the gap value for the uniform system at the 
same U. This is consistent with the local quantum criti- 
cality view [ill ll l| for which we present new compelling 
evidence in Fig. 01 Our conclusions can be tested experi- 
mentally since the dynamic structure factor S(k,oj) (and 
therefore, S(oj)) can be measured experimentally using 
Bragg spectroscopy 0,0. 
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